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From  Local  Approximation  to  a 
G1  Global  Representation 

Cedric  Gerot,  Dominique  Attali,  and  Annick  Montanvert 


Abstract.  To  represent  a complex  surface,  it  is  useful  to  describe  it  as 
a set  of  simple  parametric  primitives  such  as  quadrics.  But  if  one  wants 
to  use  few  primitives,  these  have  to  be  smoothly  blended.  To  define  this 
blending,  we  propose  to  describe  the  initial  global  surface  with  charts.  The 
blending  surfaces  result  from  a convex  combination  of  primitives  whose 
weights  are  defined  on  open  sets  of  H2  given  by  the  charts.  We  have 
established  the  properties  that  the  weight  functions  must  satisfy  to  obtain 
a G 1 representation  of  the  global  surface,  and  we  have  constructed  such 
functions. 


§1.  Introduction 

The  abundance  of  high  quality  volumetric  image  data  and  new  performant  seg- 
mentation methods  for  multidimensional  image  data  make  3-D  objects  ready 
for  analysis.  Volumetric  objects  are  basically  represented  by  a binary  voxel 
representation  or  by  a triangulation  of  the  surface.  Because  they  are  based 
on  huge  lists  of  voxels  or  surface  elements,  they  are  not  efficient  for  capturing 
global  and  local  shape  features  with  a view  to  characterizing  shape  proper- 
ties. The  spline  surfaces  can  be  very  useful,  but  become  difficult  to  use  for 
topologically  arbitrary  surfaces  modeling  because  they  require  a rectangular 
parameterization.  On  the  contrary,  any  surface  can  be  approximated  using 
quadric  surface  patches  as  in  [2],  While  they  lead  to  a good  shape  descrip- 
tion, the  quadric  patches  do  not  define  an  overall  continuous  surface.  Ideally, 
a surface  representation  for  image  analysis  should  allow  us  to  represent  con- 
tinuously any  complex  surface  with  few  parameters,  and  to  extract  shape 
properties  as  well. 

We  propose  to  represent  a surface  with  charts.  A chart  is  composed 
of  a patch  U lying  on  the  surface  and  a homeomorphism  of  U onto  a 2D- 
domain.  This  notion  has  already  been  used,  but  essentially  for  image  synthesis. 
Thus,  in  [14],  it  allows  texture  mapping  on  a triangulated  surface  which  is  too 
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complex  to  be  described  by  only  one  chart.  It  is  also  used  in  [9]  to  design  a 
surface  with  B-splines  on  any  topological  polyhedra,  and  in  [15]  to  generalize 
the  B-splines  for  constructing  surfaces  from  irregular  control  meshes  that  can 
be  embedded  in  the  plane.  Then,  Eck  et  al.  [5]  use  this  notion  to  design  a 
subdivision  mesh  from  any  triangular  mesh. 

We  use  it  for  image  analysis  because  it  allows  to  unfold  a complex  surface 
(for  instance  the  surface  of  a brain).  It  is  then  an  appropriate  tool  to  extract 
surface  features.  Before  using  it,  we  first  have  to  construct  it.  To  do  so,  we 
begin  by  representing  the  surface  by  means  of  a set  of  simple  surfaces  called 
primitives  (quadrics  for  instance)  which  approach  it  locally.  (We  currently 
study  new  processes  to  extract  primitives  from  3D  objects).  As  mentioned 
above,  the  primitives  do  not  define  a globally  continuous  surface  in  general. 
So  they  have  to  be  smoothly  blended.  This  paper  is  focused  on  a solution  to 
the  blending  problem. 

Several  different  approaches  to  surface  blending  have  been  suggested. 
Firstly  to  fill  a hole  on  a surface,  one  can  interpolate  a position  and  tangency 
conditions  network  [18,7,16,4],  or  construct  a rational  patch  to  fill  a polygonal 
hole  [10,8].  Our  blending  problem  is  not  to  fill  a hole.  But  our  approach  solve 
this  problem  too.  Secondly,  to  blend  two  surfaces,  one  can  apply  a rolling-ball 
algorithm  [1,3,6].  But  one  cannot  blend  more  than  three  surfaces  at  the  same 
time.  One  can  also  meld  isopotentials  if  the  primitives  are  implicitly  defined, 
[17,11,12].  But  this  seems  to  be  a too  restrictive  condition.  Our  approach 
differs  from  these  methods  in  that  we  blend  any  number  of  primitives  at  the 
same  time,  provided  they  can  be  parameterized. 

This  paper  begins  with  mathematical  definitions  which  are  necessary  to 
define  our  representation  with  charts.  Then  we  present  our  approach  for 
surface  blending.  Next  we  illustrate  the  different  steps  of  our  approach  with 
some  examples.  Finally,  we  conclude  with  future  work. 

§2.  Surface  Representation  with  Charts 

We  begin  with  some  mathematical  definitions  coming  from  differential  geom- 
etry [13]. 

Definition  1.  A n-dimensional  manifold  is  a topological  space  such  that  each 
point  admits  a neighborhood  homeomorphic  to  Rn . 

Definition  2.  A chart  ( U,ip ) is  composed  of  an  open  set  U of  an  n-manifold 
and  a homeomorphism  ip  of  U onto  an  open  set  of  IRn . 

Definition  3.  Two  charts  ( Ui,ipi ) and  ( Uj,ipj ) agree  with  each  other  if  their 
transition  function 

ipi:j  = ipj  o ip-1:  ipi(Ui  n Uj)  ->  ipj(Uj  n Ui) 


is  a diffeomorphism. 

Definition  4.  Such  a collection  of  maps  charting  all  of  the  manifold  is  called 
an  atlas. 
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Remark.  Every  2-dimensional  manifold  admits  an  atlas. 


Ill 


To  give  an  atlas  describing  a surface  is  to  give  a representation  of  it  with 
charts.  This  representation  possesses  two  main  advantages.  Firstly,  it  com- 
bines local  information  with  global  information  (tp~1  is  a local  parameteriza- 
tion of  the  surface).  That  means  one  can  work  locally  on  the  surface  without 
undesirable  consequences  on  the  global  surface  because  the  atlas  maintains 
this  consistency  by  definition.  Secondly,  an  atlas  allows  to  translate  problems 
given  on  any  2-dimensional  manifold  into  problems  given  on  1R2. 

§3.  Smooth  Blend 

We  have  given  the  mathematical  definition  of  the  representation  with  charts. 
We  now  discuss  how  to  use  it  to  construct  a G1  global  representation  of  a 
surface.  Our  aim  is  to  represent  a 2-dimensional  manifold  V by  means  of  a 
set  of  simple  surfaces  called  primitives  (quadrics  for  instance)  which  approach 
V locally.  To  be  more  precise,  we  assume  a family  {U,}  of  open  sets  on  V the 
union  of  which  covers  V and  such  that  each  Ui  is  approximated  by  a primitive 
Pi  in  such  a way  that  there  is  a bijection  b,  of  Ui  onto  Pi : Pi  = bi(Ui). 
The  blend  we  want  to  construct  between  the  primitives  Pi  must  be  a smooth 
surface  S which  overlaps  a closed  set  of  each  Pi,  called  pure  area  and  defined 
by  bi(Ui\ UJ5£j  Uj).  So,  to  be  able  to  construct  a smooth  blend,  the  primitives 
must  overlap  sufficiently  (see  §4.1). 

The  blend  is  defined  as  a convex  combination  of  the  primitives  P,  which 
approximate  overlapping  open  sets  Ui . The  surface  S is  defined  by  an  atlas 
and  is  a representation  of  V. 

Hypotheses. 

• We  suppose  that  Pi  are  2-dimensional  manifolds  parametrized  by  pi, 
homeomorphisms  which  are  G1  on  an  open  set  ft*  of  1R2:  P,  = pj(fij). 

• Let  Qij  be  the  open  set  of  Sfy  defined  by  fly  = p~1(bi(Ui  fl  Uj)). 

• We  suppose  that  there  exist  some  bijective  transition  functions  <fiij  : 
flij  —>  Qji,  such  that  <pkj  0 ‘Pik  = fij-  In  particular,  On  = fl,  and 
ipu  is  the  identity.  We  write  P(1N)  for  the  set  of  subsets  of  IN,  and  define 
for  all  i 

2i  :Oi  - T(N) 

m h { j £ N : mg  fly} 

• Let  the  weight  functions  a*  be  defined  on  fi,  and  satisfying  the  following: 

Property  1.  Convexity 

la)  Vi,Vm  g fl,,0  < a;(m)  < 1, 

lb)  Vi, Vm  g 0i,J2j€x i(m)  “i(fyW)  = 1- 
Property  2.  Regularity 

2a)  a,(m)  = 1 if  pi(m)  belongs  to  the  pure  area, 

2b)  ai(m)  = 0 if  m does  not  belong  to  ft,. 
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• We  define  for  all  i 

tpi  :fii  ->  R3 

^ acjfoijimflpjtipijim)). 
jeZi(m) 


Remark.  This  definition  is  consistent:  <pj(ipij(m))  — tpl{m). 

Proposition.  With  these  hypotheses  we  get: 

• If  < pi  is  bijective,  tpij  is  C°  and  Oi  is  C°,  then  S is  a 2-dimensional  manifold 

for  which  an  atlas  is  {(</:>“ 1 (Oj ) , 1 ) } - 

• If  tpij  is  C 1 and  Oj  is  C1,  then  S is  described  by  a C1-atlas  (5  is  then  a 
(^-surface). 

Property  2a  guarantees  that  <Pi(rn)  = Pi(m)  if  Pi(m)  belongs  to  the  pure 
area.  Properties  2a  and  2b  can  be  inconsistent  with  each  other  if  the  pure 
area  is  not  strictly  included  in  pi(Oj).  But  in  this  case,  another  primitive  can 
be  introduced,  which  overlaps  locally  Pj. 

On  one  hand,  our  representation  is  more  efficient  if  few  primitives  are 
used.  On  the  other  hand,  S is  closer  to  V when  more  primitives  are  used. 
So  the  appropriate  balance  must  be  found  with  regard  to  these  needs.  But, 
if  the  blend  is  not  defined  specifically  to  perform  the  approximation  of  V by 
S,  the  approximation  error  is  on  the  same  order  of  magnitude  as  that  due  to 
the  local  approximation  by  each  primitive.  This  last  property  is  due  to  the 
convexity  property  followed  by  the  weight  functions. 

To  construct  weight  functions  which  satisfy  the  convexity  properties,  we 
first  construct  functions  /3,  satisfying  the  following 
Property  3. 

3a)  Vi,Vm  € flj,0  < /?;(m)  < 1, 

3b)  Vi,  Vm  € /3j(<fiij(m))  + 0, 

3c)  Pi(m)  = 1 if  Px(m)  belongs  to  the  pure  area, 

3d)  f3i(m)  = 0 if  m does  not  belong  to  fi;. 

Then,  the  weight  functions  a,-  defined  by  the  following  expression  have 
all  the  desired  properties: 


, s _ AW 

CXiyTTLJ  . . vv  ? 

PjWji171)) 

where  ipij  is  a Cl  transition  function. 

§4.  Applications 

In  this  paper,  we  detail  the  construction  of  the  open  sets  fi,-,  the  weight 
functions  a,,  and  the  functions  ipi . Further  work  will  focus  on  the  construction 
of  domains  {/;  and  the  transition  functions 
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Fig.  1.  Functions  b(p)  and  f3(p,9). 

4.1.  Weight  functions 

As  shown  in  §3,  to  construct  satisfactory  weight  functions,  we  first  construct 
function  fa  satisfying  Property  3.  We  suppose  the  open  set  ft  is  a disc  whose 
radius  is  R.  Let  r be  the  radius  of  the  smaller  disc  having  the  same  center  as 
ft  and  including  the  set  of  points  {m  € ft  : p(m)  belongs  to  the  pure  area}. 
To  simplify  the  notations,  we  call  this  set  of  points  the  pure  area  too.  Then 
we  can  give  a cylindric  definition  of  fa  where  t is  a shape  parameter: 

/3(p,  6)  - b(p) 


with 

Kp)  = 


1 

P(P) 

Hp ) 

1 - P{R  + r - p) 

.0 


if  p <r, 
if  r < p < r + t, 
if  r + t<p<R  — t, 
if  R — t < p < R, 
if  R < p, 


where  L(x)  = Dx  + E,  P(x)  = Ax2  + Bx  + C and  D = R_^_t , E = 

A = g,  B = -2Ar,  C = 1 - Ar2  - Br. 

Fig.  2 shows  the  weight  function  a after  normalization,  in  a case  where  P 
is  combined  with  five  other  primitives. 

In  order  that  the  small  disc  whose  radius  is  r better  fits  the  pure  area, 
two  modifications  can  be  easily  implemented.  Firstly,  we  can  define  the  small 
disc  containing  the  pure  area  with  different  center  than  the  center  of  ft.  (3 
will  have  the  same  definition  but  with  R depending  on  9.  Secondly,  we  can 
use  ellipses  rather  than  discs. 

The  parameter  t,  which  belongs  to  (0,0.5),  controls  the  nonlinear  part  of 
b.  The  smaller  t,  the  smaller  this  part  is.  To  avoid  a final  surface  which  is 
visually  too  sharp,  t must  be  neither  too  small  nor  too  large. 

The  size  of  the  pure  area  also  plays  an  important  role  in  the  surface 
smoothness.  As  shown  in  Fig.  4,  if  the  pure  area  is  too  large,  then  the  transi- 
tions between  the  primitives  are  too  sharp  in  regards  with  the  resolution  of  a 
visualization  process.  On  the  contrary,  a pure  area  which  is  too  small  causes 
smooth  transitions,  but  the  shape  of  primitives  is  lost.  In  the  example  shown 


Fig.  4.  Blendings  with  too  large  (r  = 0.85H)  and  too  small  (r  = 0.07/1)  pure  areas. 


in  Fig.  3,  we  have  implemented  a case  where  the  pure  area  is  half  the  area  of 
ft.  This  balance  gives  a good  solution. 

We  have  constructed  satisfactory  weight  functions.  To  apply  our  repre- 
sentation we  must  define  the  functions  I*,  bijective  and  C1  transition  functions 
ifiij , and  check  that  ipi  is  bijective  to  be  sure  that  S is  a 6'1  surface. 
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We  first  apply  our  representation  in  the  case  where  the  surface  V to  be  repre- 
sented can  be  described  by  v(x,y)  = [x,y,  f(x,y)],  (x,y)  £ D C 1R2  : a land 
surface  in  topography  for  example. 

We  suppose  that  a set  of  open  discs  fi,  is  defined  by  any  local  approx- 
imation strategy,  and  parameterizations  pi(x,y)  = [x,y,p*(x,y)},  (x,  y)  £ 
fit  are  given  such  that  D C (Ji  and  each  pt  is  a C1  approximation  of 
{v(x,y)  : (x,y)  £ fi*  n £>}. 

It  is  easy  to  check  if  a point  (x,  y)  is  inside  a disc  f lj , and  so  to  define 
Xi(x,y).  Besides,  the  transition  functions  tfiij  are,  in  this  case,  the  identity, 
which  is  C 1 and  bijective.  Finally,  the  functions  ipi  constructed  by  convex 
combinations  of  such  pi  are  bijective.  So,  the  surface  S described  by  the  atlas 
{(<pt_1(^i),^_1)}  is  G1. 

Fig.  3 shows  an  example  of  this  first  case.  We  deal  with  six  open  discs. 
One  of  them,  fly  is  centered  on  the  origin.  The  others  are  centered  on  the 
vertices  of  a pentagon  which  encircles  Qt.  The  primitives  are  paraboloids. 
The  central  one  is  defined  by  z = k — x2  — y2,  and  the  others  by  z = (x  — 
kj)2  + {y  — k2)2  where  k*  are  constant. 

This  example  displays  the  blend  between  two  primitives  quite  similar 
locally  around  their  parts  to  blend  (a  central  and  a peripheral),  and  between 
two  dissimilar  primitives  (two  peripherals).  The  surface  is  smooth  even  if 
adjacent  primitives  are  strongly  different  from  each  other. 

4.3.  A more  general  case 

Most  of  the  surfaces  to  be  represented  cannot  be  described  by  [x,y,vz(x,  y)\. 

To  deal  with  any  surface  V,  we  require  a triangular  mesh  which  is  a first 
approximation  of  V.  To  simplify  the  notations,  we  name  this  mesh  V too.  We 
define  on  it  a set  of  domains  [/,.  Each  is  a set  of  vertices,  edges  and  faces 
of  V.  It  is  isomorphic  to  an  open  disc,  and  well  approximated  by  a primitive 
(a  plane  in  Fig.  7). 

In  this  case,  we  do  not  give  an  analytic  expression  to  the  functions  pt, 
ifij  and  the  open  set  fly,  but  they  are  defined  on  a finite  set  of  points.  They 
are  described  by  links  between  vertices  of  different  meshes  (see  Fig.  5).  For 
instance,  we  construct  a mesh  fl*  lying  on  IFt2,  using  the  bijective  harmonic 
map  presented  in  [5],  on  Ui,  see  Fig.  6.  Therefore,  each  vertex  u of  Ui  is  linked 
to  a vertex  lj  of  fly  Because  of  these  links,  2),  ipij  and  fly  can  be  defined  on 
the  vertices  of  fl,  as  follows: 

For  every  vertex  u of  V,  we  construct  l(u),  the  list  of  the  vertices 

linked  to  u.  Each  of  these  vertices  lies  on  a different  flj.  l(u)  contains 

only  one  vertex  w if  w belongs  to  a pure  area.  Let  w be  a vertex  of  an 

open  set  fly  Let  u be  the  vertex  of  Ui  (and  so  V)  linked  to  it.  For  all 

i,  if  one  of  the  vertices  of  l(u ),  w',  belongs  to  fly  then  2(u>)  includes 

j,  u belongs  to  fly  and  <Pij{oj)  = w';  or  else  u>  does  not  belong  to 
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Fig.  5.  Notation  in  the  general  case. 

The  functions  fa  are  calculated  as  in  §4.1,  and  thanks  to  Z*,  we  calculate 
on  the  vertices  of  11,. 

In  the  same  way,  we  construct  a mesh  Pi,  which  is  in  bijection  with  Ui  and 
whose  vertices  lie  on  the  primitive  which  approximate  {/;.  Thanks  to  the  links 
between  the  vertices  of  Ui  and  fi;,  we  define  the  links  between  the  vertices  of 
fi;  and  Pi . These  links  define  the  paramatrization  p,  on  the  vertices  of  fi,. 

We  then  construct  a mesh  S whose  vertices  are  calculated  by  ipi  defined 
on  the  vertices  of  Hi  as  in  §3. 

Assuming  there  exist  C 1 functions  ipij  and  pi  which  interpolate  the  values 
set  on  the  vertices  of  Hi,  and  satisfy  the  hypotheses  given  in  §3,  the  vertices 
of  S lie  on  a G 1 surface. 

Remark.  Because  we  do  not  give  an  analytic  expression  for  pi  and  ifij,  we 
have  to  store  the  meshes  Hi  and  Pi.  In  future  work,  we  will  either  have  to 
give  simple  expressions  for  those  functions,  or  decrease  the  size  of  the  meshes. 

§5.  Conclusion 

The  representation  with  charts  can  be  used  to  construct  a useful  surface  model. 
But,  before  applying  it  to  real  data,  we  still  have  two  crucial  steps:  the  defi- 
nition of  the  domains  U,  and  the  definition  of  transition  functions.  Then  we 
will  apply  it  to  image  analysis  problems  such  as  registration,  surface  feature 
extraction,  texture  mapping  or  animation.  More  precisely,  we  will  begin  with 
the  visualization  of  5 by  a mesh  hierarchy  which  offers  different  levels  of  detail. 
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Fig.  6.  An  open  set  fij. 


Fig.  7.  The  meshed  primitives  P j and  the  mesh  S. 
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